Reference Genome Research Articles

First Report of Rust Caused by Phakopsora nishidana on Creeping Fig, Ficus pumila, in South Africa.

Phakopsora nishidana has recently been reported as the causal organism of rust on edible fig, Ficus carica in South Africa (SA) (Boshoff et al. 2022). This contradicted reports by Doidge (1927, 1950) and Verwoerd (1929) who listed Cerotelium fici as the causal organism of the disease in SA. Similarities in urediniospore morphology and differing taxonomic interpretations most likely contributed to the use of both pathogen names as the causal agents for fig rust (Boshoff et al. 2022; Padamsee and McKenzie 2024). In January 2023 rust was commonly observed on nursery specimens of creeping fig, Ficus pumila, a popular evergreen and fast-growing garden plant in SA. Ficus pumila is native to southern China, Indochina, and eastern Asia but has been introduced to many countries worldwide (https://powo.science.kew.org/ accessed 20 March 2024). The F. pumila rust isolate collected in Somerset West (34°07'00.70"S, 18°52'41.75"E; Western Cape (WC) province) was identified using the sequenced 5.8S rRNA-ITS2-28S rRNA locus. When used in a BLAST analysis, this F. pumila isolate (GenBank accession number OR835538) shared the best homology with the Vermont P. nishidana isolate (34°24'33.37"S, 19°08'59.24"E) (MZ047090; e-value 0.0, 99.7% identity, 1400/1404 bp) reported on edible fig by Boshoff et al. (2022). Other P. nishidana isolates with excellent homology were KY764080 (e-value 0.0, 99.8% identity, 1012/1014 bp), KY764081 (e-value 0.0, 99% identity, 882/885 bp) and MF580676 (e-value 0.0, 99.3% identity, 981/988 bp). Accessions were used in a phylogenetic study that included four other rust samples collected from F. carica trees in Somerset West, WC (OR835534; 34°03'37.26"S, 18°51'02.18"E), Onrus River, WC (OR835535; 34°24'43.10"S, 19°09'57.99"E), Elgin, WC (OR835536; 34°14'31.37"S, 19°03'05.38"E), and Bloemfontein, Free State province (OR835537; 29°05'05.1"S 26°09'09.5"E). Reference sequences were as described in Boshoff et al. (2022). Despite limited variation, all six South African isolates, including the F. pumila isolate, grouped with the three P. nishidana reference isolates, but separate from the three P. myrtacearum isolates collected from eucalyptus trees (Maier et al. 2016). Uredinia on F. pumila leaves were mostly hypophyllous and surrounded by a halo of brown, necrotic tissue. On the upper leaf surface, the lesions appeared as small, dark leaf spots with infrequent sporulation. Urediniospores were echinulate, mostly obovoid or elliptical and their cell walls pale-yellow to yellowish-brown. The mean width and length of urediniospores were 16.7 x 21.3 µm. No telia were observed. Applying standard rust inoculation and incubation procedures (Boshoff et al. 2022), F. pumila plants inoculated with P. nishidana isolate PREM63073 produced sporulating uredinia on the abaxial surface of leaves. Likewise, inoculation of F. carica cv. Kadota leaves with the F. pumila isolate OR835538 yielded sporulating lesions containing urediniospores typical of P. nishidana, on the lower surface. Based on both the DNA sequence data and controlled infection studies, our study confirmed the host status of F. pumila for P. nishidana. The occurrence of rust on creeping fig plants in nurseries, the evergreen status of the creeper, and thus the potential to harbor the pathogen during winter, are likely to contribute to the spread of fig rust in SA.

High-quality, chromosome-level reference genomes of the viviparous Caribbean skinks Spondylurus nitidus and S. culebrae.

New World mabuyine skinks are a diverse radiation of morphologically cryptic lizards with unique reproductive biologies. Recent studies examining population-level data (morphological, ecological, and genomic) have uncovered novel biodiversity and phenotypes, including the description of dozens of new species and insights into the evolution of their highly complex placental structures. Beyond the potential for this diverse group to serve as a model for the evolution of viviparity in lizards, much of the taxonomic diversity is concentrated in regions experiencing increasing environmental instability from climate and anthropogenic change. Consequently, a better understanding of genome structure and diversity will be an important tool in the adaptive management and conservation of this group. Skinks endemic to Caribbean islands are particularly vulnerable to global change with several species already considered likely extinct and several remaining species either endangered or threatened. Combining PacBio long-read sequencing, Hi-C, and RNAseq data, here we present the first genomic resources for this group by describing new chromosome-level reference genomes for the Puerto Rican Skink Spondylurus nitidus and the Culebra Skink S. culebrae. Results indicate two high quality genomes, both ∼1.4 Gb, assembled nearly telomere-to-telomere with complete mitochondrion assembly and annotation.

A draft reference genome assembly of California Pipevine, Aristolochia californica Torr.

The California Pipevine, Aristolochia californica Torr., is the only endemic California species within the cosmopolitan birthwort family Aristolochiaceae. It occurs as an understory vine in riparian and chaparral areas and in forest edges and windrows. The geographic range of this plant species almost entirely overlaps with that of its major specialized herbivore, the California Pipevine Swallowtail Butterfly Battus philenor hirsuta. While this species pair is a useful, ecologically well-understood system to study co-evolution, until recently, genomic resources for both have been lacking. Here, we report a new, chromosome-level assembly of A. californica as part of the California Conservation Genomics Project (CCGP). Following the sequencing and assembly strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin proximity sequencing technology to produce a de novo assembled genome. Our genome assembly, the first for any species in the genus, contains 531 scaffolds spanning 661 megabase (Mb) pairs, with a contig N50 of 6.53Mb, a scaffold N50 of 42.2Mb, and BUSCO complete score of 98%. In combination with the recently published B. philenor hirsuta reference genome assembly, the A. californica reference genome assembly will be a powerful tool for studying co-evolution in a rapidly changing California landscape.

Reconstruction of 3D knee MRI using deep learning and compressed sensing: a validation study on healthy volunteers.

To investigate the potential of combining compressed sensing (CS) and artificial intelligence (AI), in particular deep learning (DL), for accelerating three-dimensional (3D) magnetic resonance imaging (MRI) sequences of the knee. Twenty healthy volunteers were examined using a 3-T scanner with a fat-saturated 3D proton density sequence with four different acceleration levels (10, 13, 15, and 17). All sequences were accelerated with CS and reconstructed using the conventional and a new DL-based algorithm (CS-AI). Subjective image quality was evaluated by two blinded readers using seven criteria on a 5-point-Likert-scale (overall impression, artifacts, delineation of the anterior cruciate ligament, posterior cruciate ligament, menisci, cartilage, and bone). Using mixed models, all CS-AI sequences were compared to the clinical standard (sense sequence with an acceleration factor of 2) and CS sequences with the same acceleration factor. 3D sequences reconstructed with CS-AI achieved significantly better values for subjective image quality compared to sequences reconstructed with CS with the same acceleration factor (p ≤ 0.001). The images reconstructed with CS-AI showed that tenfold acceleration may be feasible without significant loss of quality when compared to the reference sequence (p ≥ 0.999). For 3-T 3D-MRI of the knee, a DL-based algorithm allowed for additional acceleration of acquisition times compared to the conventional approach. This study, however, is limited by its small sample size and inclusion of only healthy volunteers, indicating the need for further research with a more diverse and larger sample. DRKS00024156. Using a DL-based algorithm, 54% faster image acquisition (178s versus 384s) for 3D-sequences may be possible for 3-T MRI of the knee. • Combination of compressed sensing and DL improved image quality and allows for significant acceleration of 3D knee MRI. •DL-based algorithm achieved better subjective image quality than conventional compressed sensing. •For 3D knee MRI at 3 T, 54% faster image acquisition may be possible.

First Report of Calosphaeria pulchella causing Canker and Twig Dieback of Peach (Prunus persica) in California, U.S.A.

California leads the United States in peach (Prunus persica L.) production, with approximately 505,000 tons produced in 2021 and valued at $378.3 million (California Agriculture Statistics Review, 2021-2022). During the spring and summer of 2023, twig and branch dieback were observed in three peach orchards (cvs. Late Ross and Starn) in San Joaquin County, California. Wood cankers and discoloration also occurred in branches, generally initiating at pruning wounds. Approximately 8 symptomatic twigs or branches per orchard were collected to proceed with the isolation of necrotic tissues on acidified potato dextrose agar (APDA). Isolations consistently yielded colonies of the fungal pathogen Calosphaeria pulchella (Pers. : Fr.) J. Schröt. (Réblová et al. 2004; Trouillas et al. 2012). Pure cultures were obtained by transferring single hyphal tips onto new APDA Petri plates. Colonies on APDA grew dark pink to red or purple in their center, with a white margin. Conidiogenesis was phialidic, producing round conidial masses at the tip of phialides. Conidia were produced abundantly on APDA, and were hyaline, allantoid to oblong-ellipsoidal, 4 to 5.5 (7) × 1.2 to 2.3 μm (n = 60). Two representative isolates (SJC-62 and SJC-64) were selected for genomic DNA extraction and sequencing of the internal transcribed spacer region (ITS) using ITS5/ITS4 universal primers and the beta-tubulin (TUB2) gene region using primers Bt2a and Bt2b. Consensus sequences of the two genes for the two isolates (ITS: PP063990, PP063991; TUB2: PP068303, PP068304) were compared to reference sequences (Réblová et al. 2015; Trouillas et al. 2012) using BLAST analysis. The ITS sequences of SJC-62 and SJC-64 were 99.8 and 99.5% identical to that of C. pulchella ex-type strain CBS 115999 (NR145357) and reference strain SS07 (HM237297); the TUB2 sequences were at least 98.5% identical to that of C. pulchella CBS 115999 (KT716476). Pathogenicity tests were conducted in 2- to 3-year-old healthy branches on 7-year-old peach trees, cvs. Loadel, Late Ross and Starn using the two fungal isolates and a control treatment (1 branch per treatment and 3 branches per tree) on each of 8-tree replicates. Branches were inoculated in June 2023 following wounding with a 5 mm cork borer to remove the bark and placing an agar plug from the margin of 10-day-old colonies on APDA directly into the fresh wound. Sterile agar plugs were used as controls. Inoculation sites were covered with petroleum jelly and wrapped with Parafilm to retain moisture. The experiment was completed twice. After four months, cankers and vascular discolorations developed around the inoculation sites. Length of vascular discoloration in inoculated branches averaged 72, 75, and 79 mm, for the Loadel, Starn, and Late Ross cvs., respectively. Calosphaeria pulchella was re-isolated from inoculated branches at 80 to 100% recovery rate, thus fulfilling Koch's postulates. The average length of vascular discoloration in the control was 13.5 mm and no fungi were recovered from control branches. Calosphaeria canker caused by C. pulchella is a global disease of sweet cherry. Recently, it was reported to cause cankers in peach trees in Chile (Grinbergs et al. 2023). To our knowledge, this is the first report of C. pulchella causing cankers and twig dieback of peach trees in the United States. These findings improve our knowledge of the etiology of canker diseases affecting peach trees and is critical for the development of effective disease management strategies.

VARista: a free web platform for streamlined whole-genome variant analysis across T2T, hg38, and hg19.

With the increasing importance of genomic data in understanding genetic diseases, there is an essential need for efficient and user-friendly tools that simplify variant analysis. Although multiple tools exist, many present barriers such as steep learning curves, limited reference genome compatibility, or costs. We developed VARista, a free web-based tool, to address these challenges and provide a streamlined solution for researchers, particularly those focusing on rare monogenic diseases. VARista offers a user-centric interface that eliminates much of the technical complexity typically associated with variant analysis. The tool directly supports VCF files generated using reference genomes hg19, hg38, and the emerging T2T, with seamless remapping capabilities between them. Features such as gene summaries and links, tissue and cell-specific gene expression data for both adults and fetuses, as well as automated PCR design and integration with tools such as SpliceAI and AlphaMissense, enable users to focus on the biology and the case itself. As we demonstrate, VARista proved effective in narrowing down potential disease-causing variants, prioritizing them effectively, and providing meaningful biological context, facilitating rapid decision-making. VARista stands out as a freely available and comprehensive tool that consolidates various aspects of variant analysis into a single platform that embraces the forefront of genomic advancements. Its design inherently supports a shift in focus from technicalities to critical thinking, thereby promoting better-informed decisions in genetic disease research. Given its unique capabilities and user-centric design, VARista has the potential to become an essential asset for the genomic research community. https://VARista.link.

Cotyledonary somatic embryo is one kind of intermediate material similar to callus in the process of in vitro tissue culture from Rosa hybrida ‘John F. Kennedy’

BackgroundRose is recognized as an important ornamental plant worldwide, and it is also one of the most widely used flowers in gardens. At present, the improvement of rose traits is still difficult and uncertain, and molecular breeding can provide new ideas for the improvement of modern rose varieties. Somatic embryos are quite good receptors for genetic transformation. However, little is known about the molecular mechanisms underlying during the regeneration process of rose somatic embryos. To elucidate the molecular regulation mechanism of somatic embryo plantlet regeneration, the relationship between the differences in traits of the two different regenerated materials and the significantly differentially expressed genes (DEGs) related to phytohormone pathways in the process of regeneration were be investigated.ResultsThese representative two regenerated samples from single-piece cotyledonary somatic embryo (SPC) culture of Rosa hybrida ‘John F. Kennedy’, were harvested for transcriptome analysis, with the SPC explants at the initial culture (Day 0) as the control. The differentially expressed genes (DEGs) in the materials from two different types for regeneration approach (SBF type: the regeneration approach type of single bud formed from SPC explants; MBF type: the regeneration approach type of multiple buds formed from SPC explants) were be screened by means of the transcriptome sequencing technology. In this study, a total of about 396.24 million clean reads were obtained, of which 78.95–82.92% were localized to the reference genome, compared with the initial material (CK sample), there were 5594 specific genes in the material of SBF type and 6142 specific genes in the MBF type. The DEGs from the SBF type material were mainly concentrated in the biological processes of GO terms such as phytohormones, substance transport, cell differentiation, and redox reaction. The KEGG enrichment analysis revealed these DEGs were more active in ubiquinone and other terpenoid-quinone biosynthesis, fatty acid elongation, steroid biosynthesis, and glycosphingolipid biosynthesis-globo and isoglobo series. In contrast, the DEGs induced by the MBF type material were mainly associated with the biological processes such as phytohormones, phosphorylation, photosynthesis and signal transduction. According to KEGG analysis, these DEGs of MBF type were significantly enriched in the porphyrin and chlorophyll metabolism, brassinosteroid biosynthesis, carotenoid biosynthesis, and peroxisome. Furthermore, the results from the phytohormone pathways analysis showed that the auxin-responsive factor SAUR and the cell wall modifying enzyme gene XTH were upregulated for expression but the protein phosphatase gene PP2C was downregulated for expression in SBF type; the higher expression of the ethylene receptor ETR, the ethylene transduction genes EBF1/2, the transcription factor EIN3, and the ethylene-responsive transcription factor ERF1/2 were induced by MBF type.ConclusionsAccording to the GO and KEGG analysis, it indicated the DEGs between two different regenerated materials from somatic embryos were significantly different which might be causing morphological differences. That was somatic embryos from Rosa hybrida ‘John F. Kennedy’ could regenerate plantlet via both classic somatic embryogenesis (seed-like germination) and organogenesis, cotyledonary somatic embryos should be considered as one kind of intermediate materials similiar to callus, rather than the indicator materials for somatic embryogenesis.

Three de novo assembled wild cacao genomes from the Upper Amazon

Theobroma cacao, the chocolate tree, is indigenous to the Amazon basin, the greatest biodiversity hotspot on earth. Recent advancement in plant genomics highlights the importance of de novo sequencing of multiple reference genomes to capture the genome diversity present in different cacao populations. In this study, three high-quality chromosome-level genomes of wild cacao were constructed, de novo assembled with HiFi long reads sequencing, and scaffolded using a reference-free strategy. These genomes represent the three most important genetic clusters of cacao trees from the Upper Amazon region. The three wild cacao genomes were compared with two reference genomes of domesticated cacao. The five cacao genetic clusters were inferred to have diverged in the early and middle Pleistocene period, approximately 1.83–0.69 million years ago. The results shown here serve as an example of understanding how the Amazonian biodiversity was developed. The three wild cacao genomes provide valuable resources for studying genetic diversity and advancing genetic improvement of this species.

A complex receptor locus confers responsiveness to necrosis and ethylene-inducing like peptides in Brassica napus.

Brassica crops are susceptible to diseases which can be mitigated by breeding for resistance. MAMPs (microbe-associated molecular patterns) are conserved molecules of pathogens that elicit host defences known as pattern-triggered immunity (PTI). Necrosis and Ethylene-inducing peptide 1-like proteins (NLPs) are MAMPs found in a wide range of phytopathogens. We studied the response to BcNEP2, a representative NLP from Botrytis cinerea, and showed that it contributes to disease resistance in Brassica napus. To map regions conferring NLP response, we used the production of reactive oxygen species (ROS) induced during PTI across a population of diverse B. napus accessions for associative transcriptomics (AT), and bulk segregant analysis (BSA) on DNA pools created from a cross of NLP-responsive and non-responsive lines. In silico mapping with AT identified two peaks for NLP responsiveness on chromosomes A04 and C05 whereas the BSA identified one peak on A04. BSA delimited the region for NLP-responsiveness to 3 Mbp, containing ~245 genes on the Darmor-bzh reference genome and four co-segregating KASP markers were identified. The same pipeline with the ZS11 genome confirmed the highest-associated region on chromosome A04. Comparative BLAST analysis revealed unannotated clusters of receptor-like protein (RLP) homologues on ZS11 chromosome A04. However, no specific RLP homologue conferring NLP response could be identified. Our results also suggest that BR-SIGNALLING KINASE1 may be involved with modulating the NLP response. Overall, we demonstrate that responsiveness to NLP contributes to disease resistance in B. napus and define the associated genomic location. These results can have practical application in crop improvement.

Unveiling the hidden diversity and functional role of Chloroflexota in full-scale wastewater treatment plants through genome-centric analyses

Abstract The phylum Chloroflexota has been found to exhibit high abundance in the microbial communities from wastewater treatment plants (WWTPs) in both aerobic and anaerobic systems. However, its metabolic role has not been fully explored due to the lack of cultured isolates. To address this gap, we use publicly available metagenome datasets from both activated sludge (AS) and methanogenic (MET) full-scale wastewater treatment reactors to assembled genomes. Using this strategy, 264 dereplicated, medium- and high-quality metagenome-assembled genomes (MAGs) classified within Chloroflexota were obtained. Taxonomic classification revealed that AS and MET reactors harbored distinct Chloroflexota families. Nonetheless, the majority of the annotated MAGs (166 MAGs with >85% completeness and < 5% contamination) shared most of the metabolic potential features, including the ability to degrade simple sugars and complex polysaccharides, fatty acids and amino acids, as well as perform fermentation of different products. While Chloroflexota MAGs from MET reactors showed the potential for strict fermentation, MAGs from AS harbored the potential for facultatively aerobic metabolism. Metabolic reconstruction of Chloroflexota members from AS unveiled their versatile metabolism and suggested a primary role in hydrolysis, carbon removal and involvement in nitrogen cycling, thus establishing them as fundamental components of the ecosystem. Microbial reference genomes are essential resources for understanding the potential functional role of uncultured organisms in WWTPs. Our study provides a comprehensive genome catalog of Chloroflexota for future analyses aimed at elucidating their role in these ecosystems.

<i>In silico</i> search for and analysis of R gene variation in primitive cultivated potato species

Pathogen recognition receptors encoded by R genes play a key role in plant protection. Nowadays, R genes are a basis for breeding many crops, including potato. Many potato R genes have been discovered and found suitable for breeding thanks to the studies of a wide variety of wild potato species. The use of primitive cultivated potato species (PCPS) as representatives of the primary gene pool can also be promising in this respect. PCPS are the closest to the early domesticated forms of potato; therefore, their investigation could help understand the evolution of R genes. The present study was aimed at identifying and analyzing R genes in PCPS listed in the open database of NCBI and Solomics DB. In total, the study involved 27 accessions belonging to three species: Solanum phureja Juz. & Bukasov, S. stenotomum Juz. & Bukasov and S. goniocalyx Juz. & Bukasov Materials for the analysis were the sequencing data for the said three species from the PRJNA394943 and PRJCA006011 projects. An in silico search was carried out for sequences homologous to 26 R genes identified in potato species differing in phylogenetic distance from PCPS, namely nightshade (S. americanum), North- (S. bulbocastanum, S. demissum) and South-American (S. venturii, S. berthaultii) wild potato species, as well as the cultivated potato species S. tuberosum and S. andigenum. Homologs of all investigated protein-coding sequences were discovered in PCPS with a relatively high degree of similarity (85–100 %). Homologs of the Rpi-R3b, Rpi-amr3 and Rpi-ber1 genes have been identified in PCPS for the first time. An analysis of polymorphism of nucleotide and amino acid sequences has been carried out for 15 R genes. The differences in frequencies of substitutions in PCPS have been demonstrated by analysis of R genes, the reference sequences of which have been identified in different species. For all the studied NBS-LRR genes, the proportion of substituted amino acids in the LRR domain exceeds this figure for the NBS domain. The potential prospects of using PCPS as sources of resistance to Verticillium wilt have been shown.

VirusPredictor: XGBoost-based software to predict virus-related sequences in human data.

Discovering disease causative pathogens, particularly viruses without reference genomes, poses a technical challenge as they are often unidentifiable through sequence alignment. Machine learning prediction of patient high-throughput sequences unmappable to human and pathogen genomes may reveal sequences originating from uncharacterized viruses. Currently, there is a lack of software specifically designed for accurately predicting such viral sequences in human data. We developed a fast XGBoost method and software VirusPredictor leveraging an in-house viral genome database. Our two-step XGBoost models first classify each query sequence into one of three groups: infectious virus, endogenous retrovirus (ERV) or non-ERV human. The prediction accuracies increased as the sequences became longer, ie, 0.76, 0.93, and 0.98 for 150-350 (Illumina short reads), 850-950 (Sanger sequencing data), and 2,000-5,000 bp sequences, respectively. Then, sequences predicted to be from infectious viruses are further classified into one of six virus taxonomic subgroups, and the accuracies increased from 0.92 to > 0.98 when query sequences increased from 150-350 to > 850 bp. The results suggest that Illumina short reads should be de novo assembled into contigs (e.g., ∼1,000 bp or longer) before prediction whenever possible. We applied VirusPredictor to multiple real genomic and metagenomic datasets and obtained high accuracies. VirusPredictor, a user-friendly open-source Python software, is useful for predicting the origins of patients' unmappable sequences. This study is the first to classify ERVs in infectious viral sequence prediction. This is also the first study combining virus sub-group predictions. www.dllab.org/software/VirusPredictor.html. Supplementary data are available at Bioinformatics online.

The genome of Citrus australasica reveals disease resistance and other species specific genes

BackgroundThe finger lime (Citrus australasica), one of six Australian endemic citrus species shows a high natural phenotypic diversity and novel characteristics. The wide variation and unique horticultural features have made this lime an attractive candidate for domestication. Currently no haplotype resolved genome is available for this species. Here we present a high quality, haplotype-resolved reference genome for this species using PacBio HiFi and Hi-C sequencing.ResultsHifiasm assembly and SALSA scaffolding resulted in a collapsed genome size of 344.2 Mb and 321.1 Mb and 323.2 Mb size for the two haplotypes. The nine pseudochromosomes of the collapsed genome had an N50 of 35.2 Mb, 99.1% genome assembly completeness and 98.9% gene annotation completeness (BUSCO). A total of 41,304 genes were predicted in the nuclear genome. Comparison with C. australis revealed that 13,661 genes in pseudochromosomes were unique in C. australasica. These were mainly involved in plant-pathogen interactions, stress response, cellular metabolic and developmental processes, and signal transduction. The two genomes showed a syntenic arrangement at the chromosome level with large structural rearrangements in some chromosomes. Genetic variation among five C. australasica cultivars was analysed. Genes related to defense, synthesis of volatile compounds and red/yellow coloration were identified in the genome. A major expansion of genes encoding thylakoid curvature proteins was found in the C. australasica genome.ConclusionsThe genome of C. australasica present in this study is of high quality and contiguity. This genome helps deepen our understanding of citrus evolution and reveals disease resistance and quality related genes with potential to accelerate the genetic improvement of citrus.

Allopolyploid origin and diversification of the Hawaiian endemic mints

Island systems provide important contexts for studying processes underlying lineage migration, species diversification, and organismal extinction. The Hawaiian endemic mints (Lamiaceae family) are the second largest plant radiation on the isolated Hawaiian Islands. We generated a chromosome-scale reference genome for one Hawaiian species, Stenogyne calaminthoides, and resequenced 45 relatives, representing 34 species, to uncover the continental origins of this group and their subsequent diversification. We further resequenced 109 individuals of two Stenogyne species, and their purported hybrids, found high on the Mauna Kea volcano on the island of Hawai’i. The three distinct Hawaiian genera, Haplostachys, Phyllostegia, and Stenogyne, are nested inside a fourth genus, Stachys. We uncovered four independent polyploidy events within Stachys, including one allopolyploidy event underlying the Hawaiian mints and their direct western North American ancestors. While the Hawaiian taxa may have principally diversified by parapatry and drift in small and fragmented populations, localized admixture may have played an important role early in lineage diversification. Our genomic analyses provide a view into how organisms may have radiated on isolated island chains, settings that provided one of the principal natural laboratories for Darwin’s thinking about the evolutionary process.

First Detection of Atractylodes Mild Mottle Virus in Atractylodes lancea (Thunb.) DC. in Hubei Province of China.

The cultivated aromatic medicinal herb Atractylodes lancea (Thunb.) DC. is widely used in the pharmaceuticals, nutraceuticals, and cosmetics industries (Na-Bangchang et al. 2014; Zhan et al. 2023). Huanggang in Hubei Province is a major production area for A. lancea (Huang et al. 2022; Wang et al. 2023). In April 2023, more than two-thirds of the surveyed plant leaves in this region exhibited virus-like symptoms, such as curling and mosaic patterns. To identify the underlying cause, 80 symptomatic plant leaf samples were collected from four fields (20 leaves per field) in this region and pooled for virome analysis. Total RNA, including ribosomal RNA, was extracted from the pooled samples using the Plant RNA Extraction Mini Kit (Onrew Biotech, Guangdong, China), for sequencing library construction. The Illumina NovaSeq 6000 platform was used to sequence the library and generate 150 bp paired-end reads. After processing the raw data with Trimmomatic software, a total of 44,354,650 high-quality clean reads were obtained. The clean reads were aligned against ribosomal RNA using BWA software (v0.7.17) to avoid interference and eliminate corresponding sequences. After removing potential contamination, contig assembly of the clean reads was performed using Megahit software (v1.2.9). The resulting contigs were compared with the virus NT database using the BLASTn program. Sequence pairwise comparison revealed 8 contigs (574 nt to 2243 nt) with identities ranging from 81.88% to 90.77% with Atractylodes mild mottle virus (AMMV, NC_027924.1, Lim et al., 2015). Additionally, contigs mapped to Carlavirus, Pelarspovirus, and other plant viruses in our virome dataset had low coverage and pairwise identity (less than 70%), which need to be further investigated. The presence of AMMV was confirmed by aligning the clean reads to the reference sequence (NC_027924.1) using BWA and SAMtools software, resulting in a consensus sequence (8024 nt) with gaps. DNA extraction from the pooled samples was performed using the Rapid Universal Genomic DNA Extraction Kit (Simgen, Zhejiang, China). Two pairs of specific primers, 3399F (5'-AAAGAAGAACCTCCTGATACGG-3')/5924R (5'-TGAACCTGATTCTCTTGGC-3') and 1830F (5'- CTCAGGAAATCCCAATGC -3')/3640R(5'-TTTCCCAATGTTCTTCGGG-3'), were designed to amplify the complete gene sequences of polymerase and coat protein (CP), based on the consensus sequence. The PCR products with the lengths of 2521 bp and 1814 bp were cloned into the pMD18-T vector (Takara Biotech, Dalian, China) for sequencing. The BLASTn analysis showed that the polymerase and CP gene sequences shared an identity of 94.51% (1929/2041 nt) and 88.41% (1419/1605 nt) with the AMMV isolate (NC_027924.1), respectively. The sequences have been deposited in GenBank under the accession numbers OR544810 and OR544811. We collected leaves from 32 A. lancea plants (16 symptomatic and 16 asymptomatic) in the fields. RT-PCR was conducted using CPF (5'-CTGCGAATATGAAAGTGC-3') and CPR (5'-GGTGAGCTTGTCTGTTAGG-3') primers, which were designed targeting a 527bp fragment of the CP gene (OR544811). Amplicons of the expected size (527bp) were detected in 24 plants (11 symptomatic and 13 asymptomatic), three of which were sequenced by Sanger sequencing, showing a 100% match to OR544811. These findings indicate that AMMV is prevalent in the major production area of A. lancea. Further research is needed to better characterize the potential risks of AMMV to A. lancea cultivation in China as well as other countries.

Telomere-to-telomere reference genome for Panax ginseng highlights the evolution of saponin biosynthesis

Abstract Ginseng (Panax ginseng) is a representative of Chinese traditional medicine, also used worldwide, while the triterpene saponin ginsenoside is the most important effective compound within it. Ginseng is an allotetraploid, with complex genetic background, making the study of its metabolic evolution challenging. In this study, we assembled a telomere-to-telomere ginseng reference genome, constructed of 3.45 Gb with 24 chromosomes and 77,266 protein-coding genes. Additionally, the reference genome was divided into two subgenomes, designated as subgenome A and B. Subgenome A contains a larger number of genes, whereas subgenome B has a general expression advantage, suggesting that ginseng subgenomes experienced asymmetric gene loss with biased gene expression. The two subgenomes separated approximately 6.07 million years ago, and subgenome B shows the closest relation to Panax vietnamensis var. fuscidiscus. Comparative genomics revealed an expansion of gene families associated with ginsenoside biosynthesis in both ginseng subgenomes. Furthermore, both tandem duplications and proximal duplications play crucial roles in ginsenoside biosynthesis. We also screened functional genes identified in previous research and found that some of these genes located in colinear regions between subgenomes have divergence functions, revealing an unbalanced evolution in both subgenomes and the saponin biosynthesis pathway in ginseng. Our work provides important resources for future genetic studies and breeding programs of ginseng, as well as the biosynthesis of ginsenosides.

High-quality faba bean reference transcripts generated using PacBio and Illumina RNA-seq data

The genome of faba bean was first published in 2023. To promote future molecular breeding studies, we improved the quality of the faba genome based on high-density genetic maps and the Illumina and Pacbio RNA-seq datasets. Two high-density genetic maps were used to conduct the scaffold ordering and orientation of faba bean, culminating in an increased length (i.e., 14.28 Mbp) of chromosomes and a decrease in the number of scaffolds by 45. In gene model mining and optimisation, the PacBio and Illumina RNA-seq datasets from 37 samples allowed for the identification and correction 121,606 transcripts, and the data facilitated a prediction of 15,640 alternative splicing events, 2,148 lncRNAs, and 1,752 fusion transcripts, thus allowing for a clearer understanding of the gene structures underlying the faba genome. Moreover, a total of 38,850 new genes including 56,188 transcripts were identified compared with the reference genome. Finally, the genetic data of the reference genome was integrated and a comprehensive and complete faba bean transcriptome sequence of 103,267 transcripts derived from 54,753 uni-genes was formed.

Genome-Wide Association Study of Body Conformation Traits in Tashi Goats (Capra hircus)

Identifying genetic markers of economically valuable traits has practical benefits for the meat goat industry. To better understand the genomic variations influencing body conformation traits, a genome-wide association study was performed on Tashi goats, an indigenous Chinese goat breed. A total of 155 Tashi goats were phenotyped for eight body conformation traits: body height, body length, chest depth, chest width, chest girth, rump width, rump height, and cannon bone circumference. Then, 100 Tashi goats were randomly selected for whole-genome sequencing and genotyped. We obtained 1676.4 Gb of raw data with an average sequencing depth of 6.2X. Clean reads were aligned to the ARS1.2 reference genome, and 11,257,923 single nucleotide polymorphisms (SNPs) were identified. The structure analysis showed that these Tashi goats were almost not genetically related. The 109, 20, 52, 14, 62, 51, 70, and 7 SNPs were significantly associated with body height, body length, chest depth, chest width, chest girth, rump width, rump height, and cannon bone circumference. Within the ±500 kb region of significant SNPs, 183 genes were annotated. The most significantly enriched KEGG pathway was “olfactory transduction”, and the most significantly enriched gene ontology (GO) terms were “cellular process”, “cellular anatomical entity”, and “molecular transducer activity”. Interestingly, we found several SNPs on chromosomes 10 and 11 that have been identified multiple times for all eight body conformation traits located in two fragments (114 kb and 1.03 Mb). In chr.10:25988403-26102739, the six SNPs were tightly linked, the TACTAG genotype was the highest at 91.8%, and the FNTB (Farnesyltransferase, CAAX Box Beta) and CHURC1 (Churchill Domain Containing 1) genes were located. In chr.11:88216493-89250659, ten SNPs were identified with several dependent linkage disequilibrium (LD) blocks, and seven related genes were annotated, but no significant SNP was located in them. Our results provide valuable biological information for improving growth performance with practical applications for genomic selection in goats.

First Report of Anthracnose Caused by Colletotrichum spaethianum on Hosta longipes in Korea.

Hosta longipes (Franch. & Sav.) Matsum. (Asparagaceae) is a perennial, herbaceous plant, native to Japan and Korea (Lee et al. 2021). In Korea, the plant is used as an edible vegetable and ornamental (Kang and Ju 2015). During 2021-2022, anthracnose symptoms were observed on leaves of H. longipes with over 70% disease incidence in Wanju-gun (35°38'47''N; 127°31'16''E) and Jangsu-gun (35°35'31''N; 127°30'03''E) in Jeollabuk-do, Korea. The disease initially appeared on old leaves, gradually spreading to young ones. The symptoms were characterized as yellow to white discoloration on the upper leaf surface with black necrotic tissue in the center of the lesion. Three H. longipes samples with anthracnose symptoms were collected. From each, a monoconidial isolate was obtained and then deposited in the Korea Agricultural Culture Collection (accession Nos. KACC 410038, 410391, and 410443). The dried specimens were housed at the herbarium of Jeonbuk National University (JBNU0129, 0137) and Korea University (KUS-F33379). Conidiomata was acervular, 65 to 80 × 56 to 70 µm in diam. Setae were dark brown, 2 to 4-septate, 63 to 161 µm long, being formed on a pale brown cushion. Conidia were hyaline, smooth-walled, aseptate, slightly curved, base truncate, 3.9 to 5.1 × 17 to 23 µm. The appressoria were solitary, olivaceous-brown, ovoid or irregularly shaped. Two-week-old colonies grown on PDA at 25 ℃ were 20-25 mm in diameter, initially white, then turned gray with age, with cottony aerial mycelium. The morphological and cultural characteristics of the fungus were consistent with those of Colletotrichum spaethianum (Allesch.) Damm, P.F. Cannon & Crous (Damm et al. 2012). To confirm morphology-based identification, the nucleotide sequences of the internal transcribed spacer (ITS) rDNA region, glyceraldehyde 3-phosphate dehydrogenase (GAPDH), actin (actA), chitin synthase (CHS1), histone (HIS3) and tubulin (TUB2) genes were determined for KACC410443, as outlined by Cannon et al. (2012) and Damm et al. (2009). The resulting sequences were submitted into GenBank (PP000829 for ITS, PP133094 for GAPDH, PP083418 for actA, PP133091 for CHS1, PP133097 for HIS3, and PP133099 for TUB2) and compared with reference sequences in GenBank using BLASTn search tool. The results showed a 100% match with C. spaethianum (MT611068), C. incanum (MN880260) and C. truncatum (EF016303) for ITS, and 100% with C. spaethianum for GAPDH (MH370513), actA (MH045677), CHS1 (MH370520), HIS3 (MH985161), and TUB (MH456884). Pathogenicity was tested by inoculating conidial suspension (1 ×104 cfu/ml) of three-week-old fungal colonies of the isolate KACC410443 onto leaves of three healthy potted plants. Prior to inoculation, leaves were deliberately wounded by pinpricking with a sterilized needle. Two wounded but non-inoculated plants served as controls. Plants were maintained in a greenhouse at 25 to 30 °C. Inoculated plants developed anthracnose symptoms after eight days, while the control plants remained symptomless. The fungus isolated from the inoculated plants was morphologically identical to that observed initially, fulfilling Koch's postulates. To our knowledge, there is no previous record of C. spaethianum on H. longipes, although C. spaethianum has been reported to infect another species, H. plantaginea (Cheon and Jeon 2016). This is the first report of this fungus on H. longipes in Korea (KSPP 2024) and globally (Farr and Rossman 2024). The anthracnose on this ornamental plant can be considered a new severe threat to planting strategies in gardens.

A phased genome of the highly heterozygous 'Texas' almond uncovers patterns of allele-specific expression linked to heterozygous structural variants

Abstract The vast majority of traditional almond varieties are self-incompatible and the level of variability of the species is very high, resulting in a highly heterozygosity genome. Therefore, information on the different haplotypes is particularly relevant to understand the genetic basis of trait variability in this species. However, although reference genomes for several almond varieties exist, none of them is phased and has genome information at the haplotype level. Here we present a phased assembly of genome of the almond cv. Texas. This new assembly has 13 % more assembled sequence than the previous version of the Texas genome and has an increased contiguity, in particular in repetitive regions such as the centromeres. Our analysis shows that the “Texas” genome has a high degree of heterozygosity, both at SNPs, short indels, and structural variants (SV) level. Many of the SVs are the result of heterozygous Transposable Element (TE) insertions, and in many cases they also contain genic sequences. In addition to the direct consequences of this genic variability on the presence/absence of genes, our results show that variants located close to genes are often associated with allele-specific gene expression (ASE), which highlights the importance of heterozygous SVs in almond.